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Abstract 

We study synchronization and noise- induced resonance phenomena in systems 
of globally coupled oscillators, each possessing finite inertia. The behavior of 
the order parameter, which measures collective synchronization of the system, 
is investigated as the noise level and the coupling strength are varied, and 
hysteretic behavior is manifested. The power spectrum of the phase velocity 
is also examined and the quality factor as well as the response function is 
obtained to reveal noise-induced resonance behavior. 
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I. INTRODUCTION 



In recent years the networks of coupled nonlinear oscillators have attracted much atten- 
tion: They serve as a prototype model for a variety of self-organizing systems in physics, 
chemistry, biology, and social sciences, and exhibit the remarkable phenomena of synchro- 
nization |T|J^. Among those the system with global coupling has been mostly studied both 
analytically and numerically, owing to analytical simplicity and some physical as well as bio- 
logical applications In such a system of globally coupled oscillators, effects of nonzero 
inertia and of noise as well as effects of periodic driving on synchronization have been exam- 
ined . Of particular interest in the presence of noise is the possible amplification of the 
response of the system, arising from the interplay between the noise and driving Such 
stochastic resonance phenomena, having various potential applications, have received much 
attention [p|-[l2|. Recently, the interesting possibility of stochastic resonance in systems with- 
out external periodic driving has been pointed out [p!^ -[T5[|. For example, noise-controlled 
resonance behavior in a periodic potential with constant driving has been discussed, and 
in the limit of low damping the inertia has been shown to play as a surrogate of external 
periodic driving . While such noise- induced resonance behavior has been investigated in 
oscillator systems with relatively few degrees of freedom, typically single oscillator systems, 
or in systems of excitable elements [|TB|, the possibility of detecting the resonance behavior 
in a system of coupled (non-excitable) oscillators has not been properly addressed. 

The purpose of this paper is to examine whether such noise-induced resonance behavior 
can appear in a coupled oscillator system with many degrees of freedom. We thus consider 
the system of globally coupled stochastic oscillators, each possessing finite inertia, subject 
to constant driving force, and investigate the behavior of the order parameter, which mea- 
sures the synchronization of the system, as the noise level and the coupling strength are 
varied. To understand the interplay of noise and driving force, giving rise to the possibility 
of noise-induced resonance behavior in the absence of periodic driving, we further consider 
the power spectrum of the phase velocity as the response to the driving force, and investi- 
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gate both the zero-frequency (dc) component and the nonzero-frequency (ac) one. The dc 
component of the power spectrum, proportional to the squared time average of the phase 
velocity, measures the probability that the system, kicked by noise, eventually escapes out 
of a potential-well minimum. Namely, it describes the inter-well transition. On the other 
hand, the ac component describes the intra-well oscillation behavior, which does not induce 
the escape out of the potential-well minimum. The order parameter is observed to decrease 
with noise, manifesting suppression of synchronization, and to display hysteretic behavior 
with the noise level as well as with the coupling strength. Suppression of synchronization 
is also reflected by growth of the dc component of the power spectrum, corresponding to 
the dispersion of the mean phase velocity, with noise. On the other hand, it is found that 
the generalized susceptibility, related to the power spectrum via the fluctuation-dissipation 
theorem, increases first as the noise grows from zero, reaches its maximum at a finite noise 
level, and eventually decreases as the noise level is increased further. Such noise-induced 
effects are also observed in the quality factor at appropriate nonzero frequencies, suggesting 
the presence of intra-well resonance. 

This paper consists of five sections: Section II introduces the system of coupled oscilla- 
tors, each possessing finite inertia, subject to random noise and constant driving force. The 
self-consistency equation for the order parameter, which measures collective synchronization 
in the system is described. In Sec. Ill, the behavior of the order parameter with the coupling 
strength and the noise level is investigated, which manifests hysteretic behavior at low noise 
levels. Section IV is devoted to the investigation of the response of the phase velocity to 
the external driving force, focusing on the interplay between noise and driving. The power 
spectrum of the phase velocity is revealed to exhibit noise-induced resonance in appropriate 
regimes. Finally, a brief summary is given in Sec. V. 
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II. SYSTEM OF COUPLED OSCILLATORS 



We begin with the set of equations of motion governing the dynamics of N coupled 
oscillators, the ith of which is described by its phase (pi {i = 1,2, N): 

■ K ^ 

fJ'4>i + 4>i + jtYI sin(0i - 4>j) =uJi + Tjiit), (1) 
j=i 

where represents the magnitude of the (rotational) inertia relative to the damping. The 
third term on the left-hand side of Eq. (|1|) denotes the global coupling with strength K/N , 
whereas the first and the second on the right-hand side represent the constant driving force 
and the random (thermal) noise, respectively. The driving force uoi is distributed over the 
whole oscillators according to the distribution giu), which is assumed to be smooth and 
symmetric about = 0. The term rii{t) represents independent white noise with zero mean 
and correlations {rii{t)r]j{t')) = 2T6ij6{t—t'), where the noise level T(> 0) plays the role of 
the "effective temperature" of the system. The set of equations of motion in Eq. (|l]) describes 
a superconducting wire network and may also be regarded as the mean-field version of 



an array of resistively and capacitively shunted junctions, which serves as a common model 



for describing the dynamics of superconducting arrays |T8[. In these cases, the constant 
driving force uJi corresponds to the direct current bias. 

Collective behavior of such an A^-oscillator system is conveniently described by the com- 
plex order parameter 

1 ^ 

vl> ^ - 5] e*^^ = Ae^^ (2) 

where non-vanishing magnitude (A ^ 0) indicates the emergence of collective synchroniza- 
tion and 6 gives the average phase. Note that the synchronized state corresponds to the 
superconducting state with global phase coherence in the case of superconducting networks 
or arrays The order parameter defined in Eq. (^ allows us to reduce Eq. (|l]) into a 
single decoupled equation 

+ 0i + KA sm{(j)i-6) =uji + r]i{t), (3) 



where A and 9 are to be determined by imposing self-consistency. Namely, the order param- 
eter, defined in terms of the phase via Eq. (0), in turn determines the behavior of the phase 
via Eq. (^), which depends explicitly on A and 9. We then seek the stationary solution with 
9 being constant, which is possible due to the symmetry of the distribution of Ui about zero. 
Redefining (pi— 9 as 0i and suppressing indices, we write the reduced equation of motion in 
the form 

li'(i) + (t) + K/\si\i(t) = uJ + rj{t), (4) 

which depends explicitly on the magnitude A of the order parameter. 

In the absence of noise (T = 0), the self-consistency equation for the order parameter 
reads 

A = - ^)mK^ + \^^9mK^f + y^g"mK^f + 0{K^f. (5) 

If the distribution g{uj) is given, the collective behavior of the system can thus be obtained 
by solving Eq. (^). In general the quadratic term of the order (fCA)^ is known to induce 
hysteresis in the bifurcation diagram [^. Accordingly, it has been concluded that the non- 
zero inertia tends to induce hysteresis in the bifurcation diagram of the system [0. 

The self-consistency equation for the order parameter in the presence of noise, particu- 
larly at such high temperatures that K/S./T <ti 1, has also been considered, yielding 



JcKiaK - 1) 



with the coefficients given by the integrals 



du g{uj) 



2(T2 + 

duj g{uj) 



4(T2 + cu2)2 8(T2 + cj2) 

6T + /i(8T2 - C^2) + ^2y(gy2 _ ^2) - 
8(r2+Cj2)(4T2 + Cu2) 

In this case collective behavior of the system has been obtained as follows: When K < 
Kc = 1/a, only the null solution (A = 0) is possible. At i^' = K^, on the other hand, the 



null solution loses its stability and the nontrivial solution A4., together with the unphysical 
solution A_ = — A+, emerges via a pitchfork bifurcation. Subsequently, it grows in a 
continuous manner {a? / y/c){K — Kc)^^^ as K is increased beyond Kc 

III. PHASE SYNCHRONIZATION 

In this section we present in detail the behavior of the order parameter with the coupling 
strength and the noise level. We have performed extensive numerical simulations on the 
equations of motion given by Eq. (|l|) at various noise levels and coupling strengths. The 
order parameter A has been computed from the definition given by Eq. (0), and its behavior 
depending on the coupling strength and the noise level has been examined. In simulations 
Eq. ([l|) has been integrated with discrete time steps of 6t = 0.001, and for convenience, a 
semi-circle distribution of radius r = 0.5 has been chosen for g{uj). (We have also consid- 
ered other types of distribution such as Gaussian, only to find no qualitative change.) In 
computing the order parameter, Nf = 10^ time steps have been used while the data from 
the first 5 x 10"^ steps discarded at each run. Both 6t and Nt have been varied to confirm 
that the stationary state has been achieved. We have then computed the order parameter 
in the system of = 2000 oscillators, each having the inertia fi = 0.8. 

The obtained behaviors of the order parameter with the coupling strength and the noise 
level are displayed in Figs. 1 and 2. Figure 1(a) shows the behavior as the coupling strength 
K is varied with the noise level T kept fixed: Circles and squares describe the behavior of the 
order parameter as the coupling strength is increased and decreased, respectively. At zero 
noise (T = 0), 20 independent runs have been performed with different initial configurations, 
over which averages are taken. The corresponding error bars have been estimated by the 
standard deviations whereas those data points without explicit error bars have errors smaller 
than the size of the symbol. Note the hysteresis manifested at zero noise and weakening 
as the noise strength is increased from zero. These characteristic features as K is varied 
for given T agree well with the results of Ref. [0. Figure 1(b) shows that the critical 
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coupling strength Kc, beyond which synchronization sets in, increases monotonically with 
the noise level T, demonstrating the suppression of synchronization by noise. Here K^. has 
been estimated by the value of the coupling strength at which the order parameter A first 
becomes nonzero to the precision of 10~^. The circles and squares, again corresponding to 
the data for increasing and decreasing K, respectively, have been obtained from averages 
taken over 10 independent runs with different initial configurations, and the error bars 
estimated by the standard deviation. Thus, unlike the excitable system [|16|, noise- induced 



synchronization does not emerge here. Further, the hysteretic behavior, reflected by the 
difference in the critical coupling strength between the two cases, is revealed to diminish 
conspicuously as the noise level is increased. 

In Fig. 2(a) the behavior of the order parameter with the noise level T for flxed coupling 
strength K is displayed. Here circles and squares represent the data for increasing and 
decreasing T, respectively, and the typical error bars, estimated by the standard deviation 
obtained from 20 independent runs with different initial conflgurations, are shown on the 
data points at K = 0.7. Figure 2(b) displays the detailed behavior with the noise strength 
for the coupling strength = 0.7, again manifesting the hysteresis. Note that the hysteresis 
is most conspicuous for the coupling strength around this value, decreasing as the coupling 
strength is increased. At weak coupling strengths, the system is not synchronized (A = 0), 
giving no hysteresis. It is thus concluded that the system exhibits quite generally hysteretic 
behavior as either the coupling strength or the noise level is varied, which has its origin in 
the non- vanishing inertia. 

IV. NOISE-INDUCED RESONANCE 

In this section we examine the phase velocity and its power spectrum, and investigate 
the possibility of the noise-induced resonance. The power spectrum of the phase velocity <pi 
is given by 
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1 ^ 



sif) = j^Y.m)r^mf)n), (7) 



i=l 



where Vi{f) = J dte^'^^^*(j)i is the Fourier component of the phase velocity at frequency / 
and the average over different noise reahzations is also to be taken. Thus ((■ ■ ■)) stands for 
the average over the noise realizations as well as over the whole oscillators. It is related to 



the response function of the system via the fluctuation-dissipation theorem [IS 



5(/)=2TRex(/), (8) 

where Re denotes the real part and the generalized susceptibility xU) is defined to be the 
Fourier transform of the appropriate linear-response function. When the system is disturbed 
by (time-dependent) external driving, the resulting change in the average phase velocity 
takes the form 

mU))) = xUWUl (9) 

where 51{f) is the Fourier component of the (uniform) external driving at frequency /. In 
particular, the dc component of the power spectrum, describing the dc response, reads 



1 Af 1 ^ r /• 

w^o)^-Ei*,(o)r^-i:/ 
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oc {m )), (10) 

1=1 ^' i=l ^ 

where (■ ■ ■) denotes the time average. In the case of a superconducting wire network or 
array, the phase velocity can be identified with the voltage via the Josephson relation, and 
the system is driven appropriately by a (time-dependent) external current. Accordingly, 
Eq. (^ connects the generalized resistance with the voltage power spectrum ||19|| . 

To investigate the dc component of the power spectrum, we begin with Eq. (^ and 
consider two types of the solution, depending on the coupling strength: In the limit of weak 
coupling strength, each oscillator in the system favors to oscillate with its own frequency 
and the system is not synchronized, yielding A ^ 0. The solution of Eq. (^ is then given 
by 

^i = u, + (t;,o - u,)e-"^' + - e-"^ f dt' e''/''r],{t'), (11) 



fx Jo 



where Vio = 4>i{t = 0) is the initial phase velocity. Taking the time average of Eq. (|TT]) in 
the stationary state {t oo), we obtain the mean phase velocity or the frequency of the ith 
oscillator, 

i^i) = uj^. (12) 
We now take the average over the oscillators; this reduces Eq. (|T2D to 

(((<^^))) = ivE('^^) = ivE^^ = o, (13) 

i=l i=l 

where the symmetry of the distribution g{uj) about = in the thermodynamic limit 
(A^ — > oo) has been used. 

On the other hand, the average of the square of the oscillator frequency, corresponding 
to the dc component of the power spectrum, does not vanish: 

• 2 1^.2 

)) = TfE('^^) =((^^'))> (14) 
i=i 

where (((^j^)) corresponds to the variance of the distribution of cjj, i.e., {{oj.-^)) = J doj g{uj) uP'. 
For example, in the simple case of the delta-function distribution gioS) = {1/2)[6{uj—ujq) + 
6{u!+llJq)], the variance is given by {{uJi^)) = uo"^, whereas for the semi-circle distribution with 
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radius r we have {{uJi^)) = Note that in this weak-coupling limit )) as well as 

{{{(pi))) does not depend on the noise level T, indicating the absence of the noise-induced 
effects. 

In the limit of strong coupling strength, the oscillators tend to oscillate in a coherent 
manner, displaying synchronization (A ^ 1). Since the order parameter A in Eq. (^ depends 
explicitly on the noise, decreasing with the noise level T, it is expected that unlike in the 

2 

weak-coupling limit {{{4>i) )) varies with the noise level. When the noise level is sufficiently 
low (T ^ 0) in this strong-coupling limit, the system is fully synchronized and described by 
the stationary solution 

0. = sin-(^), (15) 
9 



z 

which yields = and )) = 0. At high noise levels (T oo), on the other 

2 

hand, the system is not synchronized (A ^ 0), and we obtain )) = {{ooi'^)), similarly 

2 

to the case of the weak-coupling limit. Accordingly, in the strong-coupling limit, {{{4>i) )) is 

2 

expected to behave with noise as follows: At low noise levels, {{{4>i) )) increases from zero 
with the noise. As the noise level is raised further, it saturates eventually toward its asymp- 

2 

totic value, {{uJi^)). Note that {{{(pi) )) just corresponds to the dispersion or mean-square 
displacement of the oscillator frequencies over the system since {{{(pi))) = 0. Its monotonic 
growth thus indicates suppression of frequency synchronization, which accompanies that of 
phase synchronization measured by the order parameter in Sec. III. The dc susceptibility, 

2 

given by Xo = R-ex(/=0) = S{f=0)/2T oc {{{(pi) ))/T, then grows as the noise level T is 
increased from zero and diminishes with T at high noise levels; in between it is expected 
to reach its maximum. Therefore in contrast with synchronization, which is suppressed by 
noise, response of the phase velocity to the (uniform) external driving can be enhanced by 
adding an appropriate amount of noise. 

To confirm the analytical argument presented above, we have performed numerical simu- 
lations on the set of equations of motion in Eq. (0). For convenience, we have considered the 
semi-circle distribution for g{uj), and integrates Eq. (|T|) with discrete time steps of 6t = 0.01. 
In computing the phase velocity, Nt = 10^ time steps have been used at each run, with 
the data from the first 5 x lO'' steps discarded. We have again varied both 6t and Nt to 
verify that the stationary state has been achieved, and performed 10 independent runs with 
different initial configurations, over which averages have been taken. In this manner we have 

2 

computed {{{(pi) )) in the system of oscillators, for A^ up to 4096, and confirmed that there 
are no appreciable finite-size effects for A^ ^ 1000. 

Figure 3 presents the obtained behavior of the dc susceptibility or the dc component of the 
noise-divided power spectrum with the noise level T in the system of A^ = 2000 oscillators, 
each having the inertia fi = 0.8. The semi-circle distribution of radius r = 0.5 has been 
chosen for g{uj) and the coupling strength K = 3 adopted. In particular we have considered 
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both cases of increasing and decreasing the noise level, only to obtain the same results within 
error bars. The behavior shown in Fig. 3 demonstrates that noise helps the system escape 
from the potential well, enhancing the response of the phase velocity to external driving. It 
is of interest to note that Tm(~ 1-4), at which the response becomes its maximum, is almost 
the same as the critical noise strength below which synchronization sets in [Tc 1.4 for 

= 3 as shown in Fig. 2(a)]. The height of the effective potential barrier of the system 
described by Eq. @) is given by KA. Since the order parameter A decreases with the noise 
level T, the barrier height also becomes lower with T, helping the escape from the potential 
well and enhancing the response to the external driving. Eventually, at Tc the potential 
barrier vanishes and the response reaches the maximum. It is thus concluded that noise not 
only hinders synchronization, making the critical coupling strength Kc larger [see Fig. 1(b)], 
but also enhances the response of the phase velocity to the external driving force. 

We now investigate the ac components of the power spectrum, i.e., the power spectrum 
at nonzero frequencies, which gives the possibility of noise-induced intra-well resonance. For 
this purpose, we have also performed numerical simulations on the equations of motion, 
using the same parameter values, and compute the power spectrum of the phase velocity 
through the use of the fast Fourier transform. The obtained power spectrum as a function 
of the frequency / is shown in Fig. 4. At each noise level, averages have been taken over 
10 independent runs with different initial configurations, to obtain the data represented by 
such symbols as filled circles, empty circles, filled squares, etc., and the error bars have 
been estimated by the standard deviation. Note that in the absence of noise (T = 0), no 
peak appears at any finite frequencies, which is natural in the system without periodic (ac) 
driving. When small noise comes into the system, however, a peak develops at a nonzero 
frequency (/ ~ 0.3 in our simulation results) and grows up with the noise, suggesting the 
activation of intra-well oscillation by noise. As the noise level is raised, the amplitude of 
such noise-induced intra-well oscillation is expected to grow, lowering its frequency. Indeed 
the frequency at which the peak appears in Fig. 4 shifts toward lower values, demonstrating 
the noise- induced frequency shift. It eventually approaches zero frequency; this describes 
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the system kicked by noise in a potential-well minimum and escaping from the minimum. 
Namely, the intra-well oscillation induced by noise turns into the inter-well transition. To 
disclose the noise-induced effects in such intra-well motion, we have also computed the 
generalized susceptibility at several frequencies versus the noise level, where noise-induced 
enhancement in the response can again be observed. In particular, at finite frequencies, it 
is convenient to characterize such noise-induced effects by the appropriate quality factor 

Q = Smax{5 f / fmaxY^ ■, (16) 

where Smax is the peak height of the power spectrum, fmax is the corresponding frequency, 
and 5f is the half-width of the peak. Thus the quality factor Q, given by the ratio of the 



peak height to the relative width, measures the degree of the coherent motion [|T3[. We 
have computed Q from the power spectrum obtained from 10 independent runs, taking the 
average at each noise level. The obtained behavior of the quality factor Q as a function of the 
noise level T is shown in Fig. 5, which demonstrates the presence of the intra-well resonance 
induced by noise. The value T 0.7 at which Q reaches its maximum is apparently lower 
than that for the inter-well motion in Fig. 3, indicating that intra-well resonance can be 
induced by weaker noise. 



V. SUMMARY 

We have studied the synchronization phenomena and the noise-induced motion in a sys- 
tem of globally coupled oscillators, each possessing finite inertia, subject to constant driving 
force. The detailed behavior of the order parameter depending on the coupling strength 
and the noise level has been obtained from numerical simulations, which has revealed hys- 
teresis both with the coupling and with the noise as well as suppression of synchronization 
by noise. The hysteresis with respect to the coupling is most conspicuous in the absence of 
noise, weakening as the noise comes into the system; that with respect to the noise appears 
large at intermediate coupling strengths, diminishing with the coupling strength. We have 
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also considered the power spectrum of the phase velocity, as the response of the system to 
the (time-dependent) external driving, and examined the possibility of the noise-induced 
resonance in the system. The dc component of the power spectrum, which corresponds to 
the dispersion of the mean oscillator frequency, has been shown to grow with noise, again 
manifesting suppression of synchronization. On the other hand, the noise-divided power 
spectrum or the generalized susceptibility, which describes the response of the phase veloc- 
ity to the external driving, has been found to display a peak at a finite noise level, revealing 
the presence of noise-induced enhancement in the response. In particular, the noise-induced 
resonance in the intra-well motion has been observed in the behavior of the quality factor 
with the noise strength. It is thus concluded that noise in the system of coupled oscilla- 
tors not only suppresses phase synchronization but also helps the system to escape from a 
potential-well minimum in the response of the phase velocity, inducing the resonance. Such 
noise-induced resonance may be manifested by a resonance peak of the voltage power spec- 
trum in the case of a superconducting wire network. Finally, we note that the major role of 
inertia is to bring about hysteresis in the response of the system. The inertia is in general 
necessary for the system to possess (finite) natural frequencies, and expected to be essential 
to the ac resonance at these finite frequencies. On the other hand, it may not be crucial 
in the dc resonance behavior of the power spectrum of the phase velocity in the system of 
coupled oscillators. Preliminary results we have obtained for the case without inertia indeed 
indicate that the peak at zero frequency persists whereas that at finite frequency disappears, 
suggesting the presence of only the inter-well motion. The detailed investigation of this and 
other effects are left for further study. 
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FIGURES 
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FIG. 1. (a) The order parameter as a function of the couphng strength K for various values of 
the noise level T. Circles and squares represent the data for increasing and decreasing the coupling 
strength, respectively, and the solid and dashed curves are merely guides to the eye. Hysteresis is 
manifested at zero noise and observed to weaken with the noise, (b) Critical coupling strength, 
beyond which synchronization sets in, is shown to increase with T. Notations are the same as 
those in (a), with the error bars estimated by the standard deviation.. It is observed that noise in 
general suppresses both synchronization and hysteresis. 
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(a) 
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FIG. 2. (a) The order parameter as a function of the noise level T for various values of the 
coupling strength K. Circles and squares represent the data for increasing and decreasing the noise 
level, respectively, and the solid and dashed curves are merely guides to the eye. (b) Behavior of 
the order parameter with the noise level at coupling strength K = 0.7, with the same notations as 
in (a). Manifested is the hysteretic behavior as the noise level is varied. 
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FIG. 3. Behavior of the dc susceptibility xo (in arbitrary units) with the noise level T, in the 
case of the semi-circle distribution with radius r = 0.5, revealing noise-enhanced response of the 
phase velocity. Error bars have been estimated by the standard deviation and the solid curve is 
merely a guide to the eye. 
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120 




FIG. 4. Power spectrum of the phase velocity (in arbitrary units) at various noise levels: T = 
(filled circles); T = 0.1 (empty circles); T = 0.3 (filled squares); T = 0.5 (empty squares); T = 0.7 
(filled triangles); T = 0.9 (empty triangles); T = 1.1 (asterisks); T = 1.3 (plus signs). 
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FIG. 5. Behavior of the quahty factor Q with the noise level T (in the log scale), exhibiting the 
noise-induced resonance. Error bars have been estimated by the standard deviation and the solid 
curve is merely a guide to the eye. 
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